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Abstract 

We perform the two loop level renormalization of quantum gravity in 2 + e dimen- 
sions. We work in the background gauge whose manifest covariance enables us to use 
the short distance expansion of the Green's functions. We explicitly show that the 
theory is renormalizable to the two loop level in our formalism. We further make a 
physical prediction for the scaling relation between the gravitational coupling constant 
and the cosmological constant which is expected to hold at the short distance fixed 
point of the renormalization group. It is found that the two loop level calculation is 
necessary to determine the scaling exponent to the leading order in e. 
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1 Introduction 



We have developed the 2 + e dimensional expansion of quantum gravity extensively based 
on the pioneering works 0, |3|, [|, |5| and the exact solution of the two dimensional quantum 
gravity ||, [7], It is well known that the conformal invariance of the theory cannot be 
maintained in general due to the conformal anomaly. However the structure of the conformal 
anomaly is well understood in two dimensional quantum gravity and the conformal anomaly 
can be cancelled by introducing the "linear dilaton" type coupling^, |lOfl . 

We have proposed to introduce a "linear dilaton" type coupling in our formulation which 



is chosen to cancel the conformal anomaly in 2 + e dimensions |TT|, [T|, [Tj|. We have formulated 
a systematic procedure to renormalize the theory to all orders based on this idea. We have 
also developed a formal proof of the renormalizability within this scheme to all orders|T5|, IE]. 



It is important to carry out an explicit calculation at the two loop level to demonstrate the 
validity of this scheme. Our previous work in this direction is |14| in which we have computed 



the two loop counter terms to the leading order in the central charge c. In this paper, we 
complete the full two loop level renormalization of the theory by determining the counter 
terms which render the effective action finite. We also renormalize the cosmological constant 
operator to the two loop level. 

We have adopted the background gauge formalism in our investigations. We decompose 
the fields into the background fields and the quantum fields. The background fields appear 
as the external fields and the quantum fields appear in the loops. The background gauge has 
the advantage to maintain the manifest covariance with respect to the background metric. 
It enables us to use a manifestly covariant calculation procedure. We determine the two 
loop level counter terms which makes the effective action finite. The effective action is a 
functional of the background fields and hence manifestly covariant. In order to achieve this 
goal, we also need to make the two point functions of the quantum fields finite. It is because 
they appear as the external fields for the one loop subdiagrams of the two loop diagrams. 
In this paper we report the result of these calculations. 

In addition to demonstrate the validity of our renormalization scheme, there is a merit 
of these laborious calculations. They enable us to make a physical prediction for the scaling 
relation which is supposed to hold at the short distance fixed point of the renormalization 
group. In our previous works, it is found that the anomalous dimension of the cosmological 
constant operator is large which almost make it to be a marginal operator at the short 



distance fixed point of the renormalization group [11, As it turns out that the two loop 
level calculation enables us to calculate the scaling dimension of the cosmological constant 
operator and it is found to be relevant. It is because the scaling dimension is 0(e) and 
the two loop level calculation is required to determine it. We are then able to make a 
physical prediction for the scaling relation between the gravitational coupling constant and 
the cosmological constant at the ultraviolet stable fixed point of the renormalization group. 
This prediction can be compared with the scaling exponent which is measured in a recent 
numerical simulation ||17|| of four dimensional quantum gravity by putting e = 2. 

The organization of this paper is follows. In section 2, we explain our calculation proce- 
dure which is manifestly covariant with respect to the background metric. We then apply 
the covariant calculation procedure to the 2 + e dimensional quantum gravity. In section 3, 
we determine the counter terms up to the two loop level in the pure Einstein gravity theory 
with no cosmological constant. In section 4, we renormalize the cosmological constant op- 
erator and calculate the scaling dimensions of the relevant operators. We make a physical 
prediction for the scaling relation at the short distance fixed point of the renormalization 
group. We conclude this paper in section 5. The necessary information on the short distance 
singularities of the products of the Green's functions is listed in the appendix A. 

Our conventions for the geometric tensors are those of t'Hooft and Veltman[l[]. The 
background metric is denoted by and the tensor indices are raised and lowered by the 
background metric. The covariant derivatives should be understood to be taken with respect 
to the background metric. 



2 Covariant Calculation 

In this section, we describe a covariant calculation method with respect to arbitrary back- 
grounds. We subsequently apply it to the two-loop renormalization in (2 + e) -dimensional 
quantum gravity. The covariant method enables us to elucidate the general structure of the 
exact Green's functions in arbitrary Riemannian manifolds. Especially, we can express the 
short distance behavior of the Green's functions in a manifestly covariant way. This feature 
provides us a great help to perform multi-loop calculations. 



2.1 Heat Kernel Methods 



Here, we illustrate the method by which we can express the short distance behavior the 
of exact propagators in a manifestly covariant way. Although it is well explained in the 



literature [^1], [22[], we recapitulate it here for selfcontainedness. 

We consider the elliptic differential operators defined on a D-dimensional Riemannian 
manifold of the form, 

A x = I V 2 + P(x), V 2 = V^V M , (2.1) 

Here = <9 M + is the covariant derivative and P denotes a generic matrix. The exact 
propagator Ga is defined as the Green's function for the operator A in a manifestly covariant 
form. 

A x G A (x,x) = -I 5 D (x,x'), (2.2) 

We note that Ga(i,x') transforms as a bi-scalar or a bi-vector,. . . according to the types 
of the fields on which the operator A acts. 5 D (x, x') is the bi-scalar 5-function on a D- 
dimensional curved space satisfying 



d D x'y/g' 5 D {x,x') f(x') = f{x). (2.3) 
To describe the short distance behavior of the exact propagator in a manifestly covariant 



way, we adopt the heat kernel method by De Witt fL9|| . The heat kernel for the operator A 
is defined with the initial condition by 

(^T-A^ Qa{x,x'-t) = 0, 

lim Ga(x,x';t) = I5 D (x,x'). (2.4) 

We can easily write down the formal solution of the heat kernel. 

Ga(x,x';t) = < x'\ e rA \x > . (2.5) 

It is known that the heat kernel has the following asymptotic expansion in the small r 
limit. 

Ga(x, V ;t) = e"^ fX(*,z') r\ (2.6) 

(47TT) 2 n=0 

Here, a(x,x') is the bi-scalar geodetic interval defined to be equal to the half of the square 
of the geodetic distance between x and x'. We shall call it the geodetic interval. We refer 
to a multi-component function of the two spacetime points, which transforms like the direct 



product of two scalars, as a bi-scalar function. In a similar way, we define bi- vector functions 
and bi-tensor ones. a(x, x') is a single-valued function for sufficiently close points x, x' . This 
is satisfied in our case since we only have to consider the short distance behavior to obtain 
ultraviolet divergences. 

The bi-scalar of the geodetic interval is found to satisfy an important differential equation 
at an endpoint x. 



1 d 
a(x,x) = -g^{x)(T^a v , = —a(x,x'). (2.7) 



v(x,x') is the bi-scalar Van Vleck-Morette determinant, defined by 

v{x,x') = g~^(x) g~%(x') det(-a^), 

<V = dJdx^^^- (2 ' 8) 

Hereafter, we often indicate the covariant differentiations of a and v at (x') x by using 
the (primed) subscripts, v satisfies the following equation which can be derived from the 
equation (|2.7|) . 

v^V 2 a + 2a^V^ = D u3. (2.9) 

We note that a and v depend only on the background metric locally, while the Seeley 
coefficients are determined by the structure of the operator A. 

The Seeley coefficients a„ (x, x') are defined inductively by the following equations. 

a" = , a£\ = I , (2.10) 

na^ + a» = «"3 A x (v^_i), (2.11) 



which are obtained by substituting the asymptotic expansion (|2.6|) into the equations 
Here we use a vertical bar to denote the coincidence limit of the two arguments: x' = x, i.e. 
flo I = a o(%, x). The equation ( p.lOj ) shows that a^(x, x') is the two point matrix which makes 



the parallel displacement of the fields at x' to the ones at x along the geodesic connecting 
them. We also note that the a n 's depend on both the metric and the operator A in a local 
fashion. 

The equations Q2.7|), (|2.9|), ( |2.10| ) and ( [2.1 1| ) play crucial roles when we obtain the coinci- 
dence limits of (the derivatives of) a, v and a n . The coincidence limits of them are essential 
for our calculations which are based on the manifestly covariant expansions of the Green's 



functions. Here, we tabulate the results of the coincidence limits of (the derivatives of) a 
and v to the order relevant to our calculation. 
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(2.12) 
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(2.13) 



Now, we can find the short distance behavior of the exact Green's function (|2.2|) in a 
manifestly covariant form through the relation which is obtained from ( |2.2|) and ( |2.5| ). 



1 A 

/•oo 

/ dr Ga(x,x'; t) . 
Jo 



(2.14) 



Corresponding to the short proper time expansion of the heat kernel ( |2.6|) , we obtain the 
short distance behavior of the Green's function as 



Ga{x,x') = Gq(x, x') ciq(x, x) + Gi(x, x') af (x, x') + Ga{x,x), 



[2.15) 



where 



G (x, x') 
Gi(x, x') 



v 2 (x,x 



v 2 [x,x 



r(-l + f) ,„ , 



1Qtc2 



{^a^x^ + ^-a^x') 



(2-16) 



Simple pole terms in e are subtracted in the definitions of Go and G\ so that they are 
finite in the limit e — > 0. These poles reflect the infrared singularities due to the large r 
behavior of the expansion (|2.6| ). Since we have subtracted the infrared singularities of the 
Green's functions here, we can safely identify the short distance singularities of the effective 
action as the poles in e in the following sections. We have also introduced the dimensional 
regularization scale fi to ensure the consistency of the dimension. We find the well-defined 



non-singular limits of (|2.16| ) as e — ► to be 



i 

V 2 



G = ~{ log(/i 2 • 2a) + log vr + 7 } , 

47T 

l 

Gi = ^~ 2a{ log(/i 2 • 2a) + log tt - 1 + 7 } . (2.17) 

l07T 

The terms which contain logarithms cause the double-pole singularities in the two-loop 
calculations. 

In the expansion ( 2.15T) , we have separated the potentially singular pieces in the short 



distance limit a — > from the non-singular ones. Only Go and G\ produce the short distance 
singularities and they depend locally on the metric. On the other hand, the remainder G is 
regular for x ~ x' even though it contains two derivatives. It depends on the global aspects 
of the manifold in general. We note that G is also constructed to be finite in the limit e — > 
as in the cases of Go and G\. 

Combining the definitions ( [2.16| ) with (|2.9| ), Gq and G\ are found to satisfy the followings: 



(-V 2 + «"W 2 ^) Go = 5 D , 

I 

(V„ - iT*V M u*) Gi = -- o-^Go, 

(-V 2 + tr^vM) Gx = G - — (2.18) 

Air 

These relations are useful for us to derive the short distance singularities of the products of 
Go's and Gi's with certain numbers of derivatives. These singular products of the functions 
Go, Gi correspond to the short distance divergences which arise after the momentum inte- 
grations in multi-loop calculations. We perform the loop integrations in a coordinate space. 
These singular products form a basis of our calculation together with the coincidence limits. 
(For derivation, see appendix^.) 

By analytic continuation in D from D < and using the relations fl2.18| ), we find that 
for Go, 

V^Gol ~ 0, 

127re" 

Similarly, the continuation from D < 2 shows that for Gi, 



V^V.GqI ~ 7^ — R/J.V, (2.19) 



Gi| ~ V^Gi| ~ 



V^dl ~ ^-g^. (2.20) 
47re 



Here we have used the coincidence limits of a and v with the certain numbers of derivatives 
as shown in ([2.12|) and fl2.13|) . 



The coincidence limits (|2.19|) , (|2.20|) of Go an d G±, together with the ones for a n 's, easily 



lead to the limits of the Green's function. 

Ga\ ~ — 77— a>o\ + G&\, 
Zire 

V m Ga| ~ V m Ga|, 

V^Ga^I ~ --^—R^a%\ + V^V^I - g^4\ + V M G A V^|. (2.21] 
izne Zne 47re 



Here the coincidence limits of the Seeley coefficients which can be obtained from ( |2.10 ) and 
( |2.11| ) are given by 



a£| = I , 
VX| = °' 

af I = P - -RI . (2.22) 
6 v ' 

TZfj, u is defined through the commutation relation of the covariant derivatives for the fields 
on which the quadratic operator A x acts: 

[V M ,V„] = n„ v . (2.23) 

As is seen from the first relation of (|2.21|) , the coincidence limit of G appears in the finite 
non-local part of the one-loop effective action. It is also clear that the knowledge of the 
coincidence limits ( |2.21|) of the Green's function suffices to obtain the one-loop divergences. 

In general, a multi-loop amplitude is the multiple integral of a multi-component function 
of some spacetime points. However the divergent part of the amplitude after the loop 
integrations is necessarily the one of a single- component function of a spacetime point. So 
the covariant Taylor expansion |20| around a resulting point is particularly useful, when we 



obtain the divergent part of the amplitude. This is because divergent parts depend only on 
the short distance behavior of the theory. The manifest covariance in the background fields 
greatly reduces the amount of the necessary calculations. The covariant Taylor expansion is 
a series which consists of the geodetic intervals and the expansion coefficients made of the 
background fields. The expansion coefficients are usually the field strength or the derivatives 
of it. 



In dimensional regularizations, there remain only logarithmic divergences. In the theories 
near two dimensions, logarithmically divergent terms are only those of the mass dimension 
two, which are made of the background fields. Therefore we only need to expand the Green's 
functions until the expansion coefficients have the mass dimension two to obtain the divergent 
parts of the amplitude. These observations motivate us to apply the covariant calculation 
method to the 2-loop analysis of (2 + e)-dimensional quantum gravity. In the remaining part 
of this section, we spell out necessary ingredients for such an investigation. 

2.2 Applications to (2 + e)-dimensional Quantum Gravity 

Our starting point is the following action for gravitational fields, which is obtained from the 
pure Einstein action by reparametrizing the conformal mode (ip) so that its kinetic term 
becomes canonical ||13||: 



Wy = §/ d D xJg~ + m/> + ebtf; 2 ) - ^c^wj- (2.24) 

Classically parameters a and b take the values a 2 = 4eb = e/2(D — 1) while a 2 receives 
the quantum correction at the one loop level. We treat the conformal mode and the rest 
of the dee rees of freedom of the metric separately and denote the latter as g pu — g PP (& )^ v 
where h pu is a tracelss symmetric tensor. R is the scalar curvature made out of g pv . We also 
couple c copies of matter fields in the conformally invariant way as: 



' matter 



g J d D xfg {-g^d^id^ - ebtfA}, (2.25) 



The subscript % runs from 1 to c. We note that the matter fields decouple from the conformal 
mode. 

For the two-loop calculations, we need to expand the action up to the quartic order in 
quantum fields. For this purpose, we exploit the following formula in our parameterization: 

R 



+ \{h 2 Y v V li h" c V v h%+ ^(h i Y u R%+ ^V.V.^r + 0(h 5 ). (2.26) 

In the background field method, the gauge fixing terms are chosen to fix the quantum 
field gauge invariance and to be covariant with respect to the background fields. In order to 
cancel the hi/j-mixing term with derivatives, we adopt a Feynmann-like gauge: 



^ J d D xJ~g i(V„A" p - ad p ^){V v h vp - ady), 
J d D x{g I^V p h^V u h up + a^V p V u h pv 



El 

G 

+-— - — rd.ipd^ - -— ——dJjd^iA. (2.27) 
2 2(£>-l) ^ 2 2(£>-l) M J V ; 

Here the last term is a finite one-loop term necessary to guarantee the general covariance 
up to the one-loop level. This follows from the result of ref.[rB| where a 2 is related to the /3 
function of the gravitational coupling constant G. 

2 1 dG 

a ~ 2(D-1)G ^ ' 

= ^-^{e- AG +...). (2.28) 

The "linear dilaton" like coupling a is chosen to makes the conformal anomaly with respect to 
the background metric vanish. The parameter A is determined to be from the one-loop 
calculation. 

Next let us consider the gauge transformation of the theory in order to obtain F.P. 
ghost terms. The local gauge symmetry of gravity theories is the invariance under the 
general coordinate transformation. The metric changes under the infinitesimal shifts of the 
coordinates as 

Sg^ = d p e p g pu + g pp d u e p + t p d p g pv . (2.29) 

This leads to the gauge transformations of g pu , ip and (fi fields as: 

2 

dg^v = d p e p g pu + g pp d u e p + e p d p g pu - —V p e p g pu , 

5 Vi = e p d pVl + {^p^W p e p . (2.30) 
So we can derive the gauge transformation of hf 1, v field from that of g pi/ - 
5h" v = V% + V„e" + e p V p h^ - ^ v V p e" 



(2.31) 



As a result, we find the Faddeev- Popov ghost action from the gauge fixing condition 
Q2.27Q and the gauge transformations ( |2.30|) , (|2.31|) to be 



4h. 



G 



12 o 12 



^2.32) 



Here, we have used the expansion (|2.28|) of a 2 and obtained a finite one-loop term as in the 



gauge-fixing terms ( 2.27 ). 

We also need the expansion of the kinetic terms of ip and matter fields up to the quartic 
order. They are easily found, from (|2.24j ) and Q2.25Q , to be 



| J d D xfg[- l -d^d^ + -d^d^ 



^2.33) 



Finally, we show the known one-loop counter terms we have determined in |13fl which 
consist of the gravitational coupling constant renormalization and the linear wave function 
renormalization of the fields. 



let. 



25 -c + 



24?re 



6ire 
1 

127T6 



d D xJg a*l)(y P V v h» v + h? v R v 



d D xJg (V^VY + rj"R„uV 



^2.34) 



There is no renormalization of the gauge fixing term ( |2.27|) which is as expected in generic 
gauge theories. 

Now, based on the above expansions in quantum fields, we will apply the heat kernel 
method to the two-loop renormalizations in (2 + e) -dimensional quantum gravity. 



i n 



First, we collect the quadratic terms in quantum fields to define exact propagators: 

(-^gravity -^matter -^g.f. -^gh.)2? 

+ ^{(l-a 2 )v 2 + 2eW?}^ 



--^(%V 2 + 2ebR5 ij )if j 



(2.35) 



where 



/"V = ^V"* + " ( 2 - 36 ) 

is the identity for the traceless symmetric tensors in a D- dimensional curved space. 
For later convenience, we rescale the fields as 

rjr,rf^VGff,y/Grf , ^ -> . (2.37) 

So we can define the exact propagators G^^^^x, x'), G^,,(x,x'), G^/(x,x'), Gij'(x, x') for 
/i^i,, ghost, ^ and y?j fields respectively, up to 0(e), which is enough to obtain divergent 
contributions. 

A AC G c B >(x,x') = -I A b> S d (x,x') , 
A^G^V^x') = -^5 D (x,x'), 

A ik G kj >(x,x') = -6 if 6 D (x,x') . (2.38) 

Here, the quadratic operators A are as follows. 



A 



AB 



G A b'(x,x') 



( I„ 0\ 2 / 2/^ ja/3j R Q 7 V 7 V *v^(l + ^aW*^ 
V 1 / V iV2(l + ^)aI pa ^R^ 2ebR 

JabV 2 + P AB (x) , 

/ G^^a^x.x') G pv ^(x,x') \ 
\ G^p> a/ (x,x') G^(x,x') J ' 



= 9^ 2 + > 
A i:j = <%V 2 + 2ebR8 ij , 

= c%V 2 + Pij . (2.39) 



We note the traceless property of the quadratic operator ^ and the exact propagator 
G pv ^ a >[x,x') for h pu field. 

g^(x) A pv , pa = g»°(x) A„ = , 



r u (x) G^ p , a ,(x,x>) = g pW (x')G^ p , a ,(x,x') = 0. (2.40) 
These constraints originate from the traceless property of h pu field. 

r u (x) h^{x) = ^ v (x')V'(z') = o, 

G f j, v> p> er >(x,x') = <h IJV (x)hp,0i(x')> . (2-41) 



From the expansion ( 2.15| ), we should also note the traceless property of the nonlocal 
part of h pv field propagator. 

g flu (x)G^ p/a/ (x,x') = ge'°\x')G^ p , al (x,x') = 0. (2.42) 

It is clear that we can obtain similar relations for h^-ip mixing parts. 

g^(x) A^ = g^(x) A^ pa = , 

g^ix) G^(x,x') = g^'(x') G^{x,x') = 0, 

r v (x)G^,(x,x') = g^'(x')G^ p , a ,(x,x') = 0. (2.43) 

Next, let us consider the Seeley coefficients. As is seen in the previous section, the 0th 
Seeley coefficient is the geodetic parallel displacement matrix. Solving the definition 
( p.lOj ), we obtain the 0th Seeley coefficients for the fields as 

1 1 1 * / ^ / A 

«0 ^ = ~lpp'lua' + -Ipa'lup' ~ -^gpv(x)g P 'cr'{x ) , 

O-O p,f,ip' = &o ip,p'a> — , 

O-O i/n/,1 = 1 , 
gh 

a if = <V ■ (2.44) 

Here g pv is the background metric, and 7^/ is the geodetic parallel displacement bi-vector 
defined by the differential equation: 

o- p v plpul = , 

lim 7 Mi/ = g pv >{x) . (2.45) 



The bi-vector 7^/ effects the parallel displacement of a contravariant vector A u at x' to the 
covariant form at x along the geodesic from x' to x. 

Again, we note that the coefficient ao for h^ u field is traceless with respect to the 

background metric. 

g^{x) a ^ v ,p'a'{x, x') = g p ' a '{x') a ^y CT '(x, x') = 0. (2.46) 

Finally, we can find the coincidence limits of various Seeley coefficients from ( |2.39| ). They 
are the most important ingredients of our two-loop calculations, together with the limits of 
a, v and the singular products of the functions Go,G\. 
field 

V^V^Oo a/3,j'S'\ — -yplnv)ap,i'5' > 

= -^{Rfiva^gpS' + RpvaS'gp-Y' + Rfj,"l3l'9oc6' + RpufiS' #0:7' ) 5 
«1 a/3,7' <S' I — -Pa/3,7' 5' — gP Ia/3,yS' ■ (2.47) 



hu U -ip fields 



^ field 



ghost field 



fields 



a l 0/3,-0' I _ Pap,-4>' ■ (2.48) 

V^VjyOO = , 

«i = Pipip' — -R ■ (2.49) 



a a/3' I — <?a/3' ) 
V^V^ao a/3' I = -(R-fMu)a/3' , 

«1 Q /3' I = -Pa/3' ~ 7-R <?a/3' • (2.50) 
O 



a — j 

V M V y a i/| = , 

01 = Pij' — -R Sijf . (2.51) 



Here, we note that V^aol = for all fields. This is easily seen from the fact that no local 
and covariant quantity has the mass dimension one. 



3 Two-loop Renormalization of Pure Einstein Gravity 

In this section, we perform the two-loop renormalization of the pure Einstein gravity theory 
(without the cosmological constant). In order to do so, we use the heat kernel method 
discussed in the previous section. In the heat kernel method, the evaluations of the non- 
local divergences of the effective action at the two loop level are exactly the same as the ones 
needed for the one-loop quantum field renormalizations. The nonlocal divergences will be 
shown to be cancelled by introducing the counter terms we have already identified at the one 
loop level except the type which is proportional to the field equation. We then shaw that 
such a non-local divergence can be canceled by a non-linear wave function renormalization 
of h^v field. 

Our results in this section hence explicitly show that the two-loop divergences of the pure 
Einstein gravity theory can be made to be local by the one-loop quantum field renormal- 
izations and the renormalization of the gravitational coupling constant. The two loop level 
divergence of the effective action is shown to be the Einstein action form and hence can be 
canceled by the renormalization of the gravitational coupling constant at the two loop level. 

3.1 Two-loop Diagrams with ip or Matter Fields' Propagators 

First, we evaluate the two-loop diagrams with if) or matter fields' propagators. In figfj, the 
wavy and solid lines denote the propagators of h^ u and ip fields, while the dashed-and-dotted 
lines denote the propagators of the matter fields. In our previous work, we evaluated the 
diagrams of figs.|I]a-|I]d, expanding the background metric around the flat one |TJ]]. Here we 
illustrate the evaluation of the amplitudes for fig.[TJa in the heat kernel method. 
The amplitudes are given by the following three types of the multiple integrals: 

dV dV' < d ll ip{x)d ll iil){x') >< d v ip{x)d v iip{x') >< h a p{x)h a ipi{x') > 

J dV dV d^G^i d ^ ■ d v G^> d „/ ■ G a p ja ip , (3.1) 
J dV dV < d^ip(x)d^ip(x) >< d u ip(x)il)(x) >< h a / 3 (x)V v ih a i^(x') > 

, (3.2) 




e 

Figure 1: Two- loop diagrams with tp or matter fields' propagators 



1 c 



dV dV' < d p ip(x)ip(x') >< d u ip(x)ip(x') >< h a/3 (x)h a //3/(x') > R p > a i{x') 
= J dV dV' d^Gff ■ d v G^i ■ G a /3 >a ip>Rpicri , (3-3) 

where dV denotes an invariant volume element dPx^fg. 

We can evaluate the divergent parts of the above three types of amplitudes for general 
subscripts. We first substitute the short distance expansion (|2.15| ) of the exact propagators 



into the equations (|3~T|)- (|0| ). The expansions of them are found to consist of the products 
of the functions Go? G\ an d the Seeley coefficients and G. In order to extract divergent 
contributions, let us look at the products of the Seeley coefficients and G. We only have 
to keep the products of the quantities with the mass dimensions less than or equal to two, 
when we make the covariant Taylor expansions of them. This is seen from a simple power 
counting argument. 

Furthermore, we have to extract the singular behavior proportional to a delta function 
from the products of the functions Go, G\. (See appendix 0.) The singular behavior with a 
delta function corresponds to the short distance singularity due to the loop integrations. In 
this way, we obtain the products of the necessary Seeley coefficients and a delta function as 
the short distance behavior of the integrands of the multiple integrals. The delta functions 
with no derivative make it possible for us to take simple coincidence limits of the Seeley 
coefficients and G. But in the case of the delta function with the derivatives, we have to pay 
a little more attention to take coincidence limits. For instance, we may use an identity as 

fV 2 5 D = V 2 (/| • 5 D ) - 2V(V„/| ■ 5 D ) + V 2 /l ■ S D . (3.4) 

In such a procedure, we can obtain the singular parts of the amplitudes (|3.1|)-( |3~3"D . 

As an example, we can write down the explicit form of the singular contributions of the 
amplitude ( |3.1|) as follows. 



d D V dPY < d^(x)d^ip(x') >< d u ip{x)d v ,ip(x') >< h aP (x)h a ' pf {x') > 

~ J d D V d D V [(5 Vrt of i/d^t p , ■ <y/ d v G Q *d j ■ G ) g f V#V W/?' 

+ (5 Vrt of 7/^0 d p , ■ j/d u G 8 a , ■ Gi) 

+ (5 D part of d pl ■ y/dvGi d a , ■ G ) 

+ (5 Vrt of 7 /<%G d „ ■ j/d^Gt d a , ■ Go) Q P v>g a n'PipipIafaa'P' 

+ (V M V,(5 D part of 7/ 'd li G > d ' ^ ■ 7/ d v G *d 



with V P VJ U — > 5 U ■ g%'9\'{y^uG a p )alpl \ + l)/2 



+ (5 part of 7 P P 9 M G d p , ■ 7 CT CT d a >) g p ' p ,g° '^G^^p \ 

+ (5 D part of 7/^0 d p , ■ 1 / d v G x d „,) 

+ (5 D part of 7/ d v G Q d p > ■ ^/ d^Gxd a >) g p v>g a ^P^G^^p^ 

- (5 D part of 7 P ^ 'd p G *d ' p , ■ G ) g p p >I a p,a>p> Vp V V G^\ 



(3.5) 



- (5 part of 7/ d u G d p < ■ G ) g P u>Iap,a>p'^ V M G^| 

Here, the 5 D part denotes the singular one of the product of the functions Go,G\, which 
is proportional to a delta function. Also, a quadratic derivative part of a delta function is 
meant by the V P V U 5 D part. 

As a result, the total divergent amplitude for fig.[I]a is found to be 
G f , D m -7 + 966 
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Here, we note that there arise only the single-pole divergences. In general, the two-loop 
quantum corrections give rise to the double-pole divergences as well as the single-pole ones. 
However, the cancellation of the double-pole divergences is necessary to ensure the consis- 
tency of the theory, such as the finiteness of the (3 functions of the renormalization group. 
In the remaining part of this section, we will also find each type of the diagrams to yield 
only the single-pole divergences. 

Let us proceed to the evaluation of the amplitudes of fig.|T]b. They are of the form of the 
square of the one-loop amplitudes. Here, we also have the three types of the expectation 
values to consider. 



d D xy g < d p ip(x)d v ip(x) >< h a p(x)h 7 s(x) > 
J d D x\Jg~ d p G 1 pif,d v \ ■ G a p a §\ , (3.7) 
J d D x^J~g~ < V p h a(3 (x)V v h 1& (x) >< ip(x)ip(x) > 
J d D x\J~g~ V ij,G a p^s^ '' v\ -Gjfnp] , (3.8) 
d D xJg < h a p{x)h lS {x) >< vfj(x)ip(x) > R pi , pa (x) 



d x\Jg G a p n s\ ■ G^l ■ R^vpa . (3-9) 

The evaluation of the above is much easier than that of figJIJa. It is straightforward to 
perform it by using the coincidence limits ( [2.21 ) of the exact propagators. 



As an example, we explicitly write down the divergent contributions of the amplitude 



d D xy g < d /J ip(x)d l/ ip(x) >< h a!3 (x)h 1 s(x) > 
^ 2 
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We have obtained the total divergent amplitude for fig.[l|b as follows. 



(3.10) 



+ ^ / dPxfg Q-R, V G%^\ - ^RG^ IV \ ~ ~V 2 Sy) • (3.11) 

We are left to consider the amplitudes for figs.|I]c,[I]d. There is little difference between 
if) and matter fields except for the term: atpR. (See the actions (|2.24j ) and (|2.25|) .) It is 



irrelevant to the evaluation of these amplitudes, since this term is linear in ip. As a result, 
we have found that the amplitudes for figs.|l]c,[l]d are given by (|3~6| ) and (|3.11|) , respectively, 
multiplied by the number of matter fields c and with G^ replaced by G w . 

The amplitudes for fig.p]e give neither double-pole nor non-local divergence. This is 
because the products of the three-point vertices are proportional to a 2 ~ e/2. Although a 2 
receives the quantum correction, it can be ignored here since it is at higher orders in G. In 
this way we have found the divergent amplitude of this diagram to be: 

G r d D xJ§ (--)r. (3.12) 
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Here, we note that this result is obtained from the double-pole divergence suppressed by 
a 2 ~ e/2. Such underlying double-pole singularities arise only in (|3.12| ) and ( |3.17|) . 

3.2 Two-loop Diagrams with Only Field Propagators 

Secondly, we consider the two-loop diagrams made only of the propagators of field. 



a 



b 



Figure 2: Two-loop diagrams made only of h^ u field propagators 

In a quite similar way to the previous subsection, we can evaluate the amplitudes for 
fig.Qa. We only have to consider three types of the expectation values. We can also evaluate 
the divergent parts of them for general subscripts. Using them, we have found the total 
divergent amplitude for fig.^a to be 

G r , n r. 15 



(to) 



6 - — 49 - — 31 \ 

+-R w »*G^o°\ + —R^G^l - ^- e R G^r\) ■ (3-13) 

As for the amplitudes of fig.||b, there are two types of the expectation values to be 
considered. We have obtained the total divergent amplitude for fig|§b as follows. 

G r ,n n ( 11 
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+ §;J d "^J (-|v 2 G^,""| + |v„V„GV1 

--^C^-^l - — iVGV'| + i^^ri) • (3-14) 

3.3 Two-loop Diagrams with Ghost Propagators 

Thirdly, we consider the two-loop diagrams with the ghost propagators. In fig.^|, the dashed 
lines denote the ghost propagators. 

We have three types of the diagrams figs.|3]a-[3]c for the ghost contribution. In a similar 
way, we can evaluate these amplitudes. The total divergent amplitude for fig.^a is 

G f , n m 11 
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Figure 3: Two-loop diagrams with ghost propagators 
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For fig.|3|b, we have found the sum of the divergent amplitudes to be 
G r .n /7 / 14^ 
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+£v 2 G§ h /l - £^^ gh ^|) . (3.16) 

The amplitude for fig.|3]c gives only the local single-pole divergence. This is also because 
the product of the three-point vertices is proportional to a 2 ~ e/2. The result is 

°' ^ d D xJ~g (~)R • (3.17) 
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Again, we note that the result is obtained from the underlying double-pole singularity sup- 
pressed by a 2 ~ e/2. 

3.4 Contribution from The h^-t/j Mixing Terms 

Here, we consider the diagrams with the mixings of and i/j fields' propagators. 

There are four types of the diagrams figs.§a44|d which give nontrivial contributions. Each 
of them has only single mixing of and ip fields' propagators. The diagrams with the 
double mixings or more are found to give no divergent contribution. 

The evaluation of the diagrams with the mixings is much easier than that of the ones 
without the mixings, since there does not exist the parallel displacement matrix clq connecting 
hfj, v field to ip field. We have obtained only the non-local divergences for each of the diagrams. 



This is because the local divergences necessarily arise from the coincidence limit ( 2.48|) of 
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Figure 4: Two-loop diagrams with the mixings of the h^ u and ifj fields' propagators 

. \ 



Figure 5: Two-loop diagrams with the insertion of one- loop finite terms 

the 1st Seeley coefficient a\ but it is a traceless matrix. As a result, we have found the 
sum of the amplitudes to be the following non-local divergences. 

^ J d D xfg ^at(V,V u G^^\ + iVG^I) . (3.18) 
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As it is expected they are just cancelled by the contribution from the counter terms ( |3.20| ). 



3.5 Contributions from The One-loop Finite or One- loop Counter 
Terms 

Finally, we evaluate the amplitudes for the diagrams with the insertion of the one-loop finite 
terms or the one-loop counter terms. 

We have the only diagram fig.[| to consider for the contribution from the one-loop finite 



a 



b 



Figure 6: Two-loop diagrams with the insertion of one- loop counter terms 



terms. The divergent part of the diagram is found from (|2.32| ) to be 
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lD n -25 + c 
a x\ I g — R . 
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In a similar fashion, we can easily evaluate the diagrams with the insertion of the one- 
loop counter terms. We have obtained the sum of the divergent amplitudes for figs.^a and 
|b from ( 12T2T| ) and ( gg§ as follows: 
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+ I(v 2 G sh /l-^ gh ""!)} . 
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(3.20) 



3.6 Total Amplitude 

After these calculations, the sum of the two loop divergences except those from the counter 
terms is evaluated by summing the contributions of ( |3.6| ), ( |3.11| ), ( |3.12j ), ( |3.13| ), ( |3.14| ), ( |3 . 1 5| ) , 
Q)> 0> (EH ^d (£JJ) : 

G f aD fz 113 + 3c 

W) 



f , n r 113 + 3c - 

/ d D xJg — R 

J V y 72e 
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Here, we note that the non-local divergences originate from the one-loop sub-divergences mul- 
tiplied by the non-local parts of the propagators with possible derivatives G, V 2 G, V^V^G. 
This is directly seen from our calculation procedure. When we evaluate a non-local diver- 
gence at the two-loop level, we replace one of the propagators of a two-loop diagram by its 
non-local part G. A non-local divergence is just a non-local part G multiplied by a local 
divergence of a sub-loop which exist in the relevant diagram. 

Therefore the non-local divergences at the two loop level exactly match the one-loop 
divergent corrections of quantum two-point functions. For example, let us consider a pair 
of the three-point vertices of quantum fields. We choose a pair of quantum fields from the 
each of the vertices as external fields. Next, we make a loop connecting the remaining two- 
point vertices. The one-loop divergent correction of a quantum two-point function is the two 
external fields multiplied by the local divergence of the loop. This is identical to a non-local 
divergence of a two-loop diagram, if we replace the two external fields by the non-local part 
of their propagator G. Conversely we have the following transformation rules, which give 
a one-loop divergence for a quantum two-point function from a non-local divergence of a 
two- loop diagram. For example, 

e e 
C ^R^G^\ — ^LR^hrhT. (3.22) 

We can therefore determine the one-loop quantum field renormalizations which are re- 
quired to cancel the nonlocal divergences of the two loop amplitudes if we evaluate the 
non-local parts of the two-loop corrections. The renormalizability of the theory implies that 
the necessary counter terms can be supplied by the renormalization of the couplings and 
the fields. On the other hand, we have already determined the counter terms (|2.34j ) which 
consist of the one-loop renormalization of the gravitational coupling constant and the linear 



wave function renormalization of the quantum fields|13 



The sum of the amplitudes ( |3.21| ) and (|3.20| ) are indeed found to be local except the term 



which is proportional to the field equation: 



G f jD [7 -79 + 3c 6 
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+ £j d D xfg - \g„R)G%r P \ ■ (3-23) 

Since the field equations appear when we take the variation of the action by the fields, 
the remaining divergence can be dealt with by a nonlinear field renormalization at the one 



loop level. Previously we missed this counter term since we considered only the linear wave 
function renormalization of the fields. We remark that there is no principle to prohibit the 
nonlinear field renormalizations in this theory since the fields are dimensionless. Therefore, 
we are led to introduce the following renormalization of field, paying attention to the 
traceless property. 

h% ={l- — ) V + Z n . L (h w h p u - lg^h pa h p ^ . (3.24) 

From the linear terms of the expansion ( [2.26|) we obtain additional counter terms, and find 
the expectation values of them to give only the non-local divergences: 



J d D xfg 2Z n ,.(iV - l -g, v R)G%^\ . (3.25) 



Therefore, the two-loop expectation value (|3.23 ) can be made to be local, only if we choose 



the non-linear renormalization constant Z a i of h pv field as 

z - = • (3 - 26) 

Finally the total divergent contribution to the effective action at the two loop level in 
pure Einstein gravity is found to be of the renormalizable form. 

G ' d o x Jj - 79 + 3e fl. (3.27) 



(Air) 

Here, we note that there arises only the single-pole divergence. The cancellation of the 
double-pole divergences is necessary for the consistency of the theory. 

4 Renormalization of the Cosmological Constant Op- 
erator 

In this section, we study the renormalization of the cosmological constant operator. As 
it will be shown, we need to introduce new counter terms to renormalize the cosmological 
constant operator. We also show that these new counter terms do not spoil the result of the 
previous section. We consider an infinitesimal perturbation of the theory by the cosmological 
constant operator by adding the following term to the action: 

A J d D x^g = A J d D xfg{l + — ^ 

~ A / d D Xy f§ exp((l - - + •••)• (4-1) 

J 2a 8cr 



This term is manifestly generally covariant and it is invariant under the gauge transformation 
( p.30| ). The quadratic part of the action with respect to ip field becomes: 

/ dDx ^>{~2G {l ~ ^w* 1 * + A ^ - + h ehk ^) ■ (4 - 2) 

We regard A to be infinitesimally small and evaluate the divergent part of the effective 
action proportional to A. Namely we consider the diagrams with the single insertion of the 
cosmological constant operator. Further insertions of the cosmological constant operator 
do not lead to the short distance divergences by the power counting. We utilize the same 
covariant calculation method which has been employed in the previous section. The relevant 
diagrams in this section are those which contain if) field propagators. The divergent part 
proportional to A at the one loop level is 

- / d D xfg{^-)k{ ± + I) . (4.3) 
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At the two loop level, we find the following divergences: 

Here we omit the square of the one loop result ( |4.3|) . In this expression, we find two types 
of the nonlocal divergences. 

The first type arises due to the one loop level divergence of ^A / d D Xyfgh llv h llv form. 
It also leads to the double pole singularity in (|4.4j). Such a divergence can be subtracted 



by the nonlinear wave function renormalization of if) field as ip° = if) + a(^-)(h flv h tJ ' u + — ). 
The last term at 0(G 2 ) is added in order to cancel the leading singularity of the vacuum 
expectation value of if). By doing so, we can eliminate the first nonlocal divergence in ( |4.4j ) 
and obtain 

/^(( ° + ^ + (i)(^)) + |VV^,|) . (4.5) 

This procedure also gives rise to the following counter term from the Einstein action: 

(-£-) / d D xfgR(h, u h^ + -) . (4.6) 
47re J v ixe 

At the one loop level, it can be regarded to be a finite counter term by expanding a 2 in terms 
of G. It leads to no new divergence of / d D x\fgR type at the two loop level. Therefore this 
new counter term is harmless to the order in the perturbation theory we are working here. 



We have also found the following contribution of the ip field kinetic term type due to 
the one loop subloops of h^ v and ghost fields which is responsible for the second nonlocal 
divergence in 
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|rr . , d u x^gd^d^ . (4.7) 

This one loop contribution is another source of the double pole singularity of the cosmological 
constant operator at the two loop level. It is because a 2 factor in the numerator of this term 
can be cancelled by the inverse power of a which is associated with if> field in the cosmological 
constant operator. We need to subtract this term from the one loop subdiagrams when we 
consider the diagrams with the single insertion of the cosmological constant operator. In 
order to do so, we have chosen to adopt the following counter term in the action: 

2ne 



d D x^(^R(l + aip + ebi) 2 ) - -d^d^ip) . (4. 



This action contains the desired counter term of the ip field kinetic term type. Although it 
contains other counter terms, they can be regarded as finite counter terms by expanding a 2 
in terms of G. It is because ip field appears with additional factors of a or e as long as they 
are concerned. The terms which involve only h^ u field are safe since they do not couple to ip 
field directory. After this procedure, the two loop level divergence proportinal to A becomes: 



The counter terms (|4.8| ) do not give rise to new divergences of the Einstein action form 
neither at the two loop level. 

After these considerations, we find that the divergence of the effective action which is 
proportinal to A is local. However it still contains the double pole in e. The double pole 
here can be removed by the following wave function renormalization. We note that there is 
a freedom of the scale transformation of the background metric. We utilize this freedom to 
eliminate the divergences of the cosmological constant type which is independent of a 2 . We 
rescale the background metric by the following factor: 

After the scaling, the remaining two loop level divergence proportinal to A has only the 
simple pole in e: 



^(§) 2 A((^)(±) +( i)(^)). (4.11: 



This procedure also leads to the rescaling of the Einstein action. 



Z* J d D x^-^R 
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Here, the exponent: | is seen from the dependence of the Einstein term and the cosmological 
constant operator on the conformal mode. 

y/gR-e-it-JtjR , y/g~ e -%+J§. (4.13) 



The rescaling of the background metric introduces a new counter term as in Q4.12| ). Due to 



the new counter term, the total two-loop divergence of the effective action which is of the 
Einstein action type becomes the following instead of the final result of the previous section 

° U° X Jl ~^- c) R. (4.14) 
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It has a simple pole in e and vanishes when c = 25 which corresponds to the critical string 
theory in the two dimensional limit. These are certainly the desirable features. 

Here we recall the counter terms at the one loop level which can be associated with the 
bare gravitational coupling constant 

li%25-c t (D-l)a 2 
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The first term is the coefficient of / dPxy/gR and the second term is the coefficient of the 
Einstein action in ( |4.8| ). The only first term contributes to the conformal anomaly at the 
one loop level and enters in the definition of a 2 at 0(G) [0. We note that the second term 
contains 0(G) term. We have to add an appropriate 0(G) term to it in order to cancel 

(EH)- 

In this way, the bare gravitational coupling constant at the two loop level is found to be 
1 /jf G / ,25 -c w G N ,5(25 -c) w G 2 



We can still remove the factor (1 + |^) from this quantity by rescaling the background metric. 
We also need to multiply the factor (1 + ^)~^ to the cosmological constant operator at the 



same time. In such a renormalization scheme, the bare gravitational coupling constant is 
found to be: 

1 n € f^ ,25 -c w G\ ,5(25 -c) w G 2 



The (3 function of the gravitational coupling constant is 



Po = eG- ( 2 \^)(^)G - (*^^-)(£-)*G (4-18) 



It agrees with our previous result of to the leading order in c. 

The remaining simple pole divergences of the cosmological constant operator which are 
inversely proportional to a 2 can be ascribed to the singular expectation value of ip 2 '- 

<^>=^(l + (7)^ ( 4 - 19 ) 
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They are of order (G/a 2 ) and they can become large around the short distance fixed point 
of the renormalization group where a 2 — > 0. At higher orders, we also find the divergences 
of order (G/a 2 ) n . Therefore we need to sum this class of contributions to all orders. It can 
be done by considering the following integral [pT| . 

J d^exp(Uog{l + - ^(1 - (\)^ 2 ) (4.20) 

This integral can be evaluated for small e by the saddle point method with the change of 
the variables: ^ = p. The saddle point of p satisfies: 

Po(l + Po)-^(l + (7)^=0 (4-21) 
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The renormalization factor is found to be 

/4, . . 167ra 2 . G . 9 \ G . , 

Z A = exp^-log(l + p Q ) - -^-(1 - I^)/ > oJexp(- — ) (4.22) 

where the last factor is necessary to cancel the one loop divergence which is independent of 
a 2 . The anomalous dimension is found to be 

7 = V-0-logZ A , 
op 

We therefore find the anomalous dimension to be 



G 8tt / „ / A Go 2 , G , 



In these considerations, a 2 is regarded as 0(G) quantity and the anomalous dimension is 
0(1) at the one loop level. At the two loop level, we find a finite correction of 0(G) to the 
leading result. 

At the short distance fixed point of the renormalization group where a 2 vanishes, we find 
the anomalous dimension of the cosmological constant operator to be 

2(1 - — ) (4.25) 

where G* = + 0(e 2 ) is the fixed point coupling constant. The scaling dimension is 
the sum of the canonical and the anomalous dimensions. The scaling dimension of the 
cosmological constant operator around the fixed point is found to be 

_ (2 + £) + 2 _£: = _ £ __il_ + 0(£ *) (4.26) 

The leading order quantity of 0(1) has been cancelled by the quantum correction and the 
scaling dimension is found to be 0(e) by the two loop level calculation. 

The scaling dimension of the inverse gravitational coupling constant around the fixed 
point is: 

4^-^ )= -4-^ +o ^ < 427 » 

Since the scaling dimension of the cosmological constant operator is determined to be 0(e) by 
the two loop calculation, we only need the one loop (3 function of the gravitational coupling 
constant here. It is because of the cancelation of the leading 0(1) term in the scaling 
dimension of the cosmological constant operator. In order to determine the scaling dimension 
up to 0(e 2 ), we need to renormalize the cosmological constant operator up to higher orders. 
We find that the both gravitational and the cosmological couplings are relevant. Although 
they are also relevant couplings classically, the scaling dimension of the cosmological constant 
operators has changed drastically from the classical value of D to that of 0(e). However 
these scaling dimensions separately cannot be physical quantities. It is because they are 
scheme dependent on how we rescale the metric. 

The physically meaningful prediction is the scaling relation between the gravitational 
coupling constant and the cosmological constant. We find the scaling relation between the 
gravitational coupling constant and the cosmological constant at the short distance fixed 
point as 

" ^r + ^ ~ ^ (4.28) 



It is very different from the classical scaling relation which holds at the weak coupling infrared 
fixed point: 

~ A* (4.29) 

These scaling exponents have been numerically measured in four dimensions by using 
a real space renormalization group technique ||17||. What has been observed is the following 
scaling relation at an ultraviolet stable fixed point: 

ft) 

N~- ~ \ K * - k\ (4.30) 

where k is the inverse gravitational coupling constant and N is the space time volume which 
is dual to the cosmo logical constant. Therefore it leads 

A^~|k*-k| (4.31) 



(3q is measured to be 6.5(1.2). If we put e = 2 in our scaling relation (|4.28| ), it gives (3q ~ 4 



On the other hand, the classical scaling relation corresponds to the exponent (3q = 2. We 
observe that our prediction is closer to the measured value than the classical exponent. 

5 Conclusions and Discussions 

In this paper, we have performed the two loop level renormalization of the quantum gravity 
in 2 + e dimensions. We have adopted the background gauge and computed the effective 
action of the background fields which is manifestly covariant. We have made the most of the 
covariance to simplify the calculation. We have shown that the short distance divergences 
are renormalizable within our scheme to the two loop level. 

In fact we have explicitly determined the counter terms which make the effective action 
finite up to this order. These counter terms are all of the expected types which can be 
supplied by the renormalization of the couplings and the wave functions. The most nontriv- 
ial counter term is the type of J d D x\fgR which gives rise to the conformal anomaly. Our 
renormalization scheme is successful to handle it by choosing the "linear dilaton" type cou- 
pling a in such a way to maintain the conformal invariance with respect to the background 
metric. The bare action can be chosen to be a manifestly covariant form simultaneously |T6| . 



We have thus demonstrated the effectiveness of our formalism by the explicit calculations 
reported here. 



We have further predicted the scaling relations and the scaling exponents between the 
gravitational and the cosmological coupling constants. Such a prediction can be compared 
with recent numerical investigations in quantum gravity. We have found that our prediction 
appears to be closer to the measured value than the classical prediction. However it also 
turns out that taking the continuum limit is not straightforward at the ultraviolet stable fixed 



point where the scaling is observed (18]]. The problem is that the phase transition appears 
to be of the first order rather than the second order. So constructing a realistic quantum 
gravity theory still is a challenging task and the exponent which has been measured in the 
numerical approach may be taken qualitatively at best. 

However we cannot rule out the possibility that the continuum limit may be taken with 
further ingenuity and the scaling behavior currently observed may qualitatively survive. 
Therefore we believe that it is worthwhile to seriously attempt to calculate the scaling 
exponents in the continuum approach in view of the progress which has been already achieved 
in the numerical approach. We hope that our results will shed light on the continuum 
quantum gravity with such an optimism. 
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A Singular Products of The Functions Go, G\ 

In the course of the calculation of the two-loop amplitudes, we need to know the singular 
behavior of the several products of the functions Go's and G\% at short distance. In general, 
the singular part is given by the bi-scalar delta functions with derivatives or those multiplied 
by covariant curvature tensors. In terms of them, we can evaluate the loop amplitudes in a 



manifestly covariant way. In this appendix, we will illustrate the derivation of them [23] . 



The functions G and G\ are effectively defined by the equations ( |2.18| ). The consistency 



with the first of ( 2.18 ) and the coincidence limits (|2.19 ) lead us to find the singular behavior 



V,V,G • G " ~ - l D } { l e)n g,J D , (A.l) 

from which we can immediately obtain, discarding a non-singular total derivative, 

V M G V,G ■ Go"" 1 ~ — To v» ■ (A.2) 

Differentiating ( |A.1| ) and ( |A.2| ), we can derive the singular behavior of the products with 
three derivatives. They are proportional to the delta functions with a derivative, since there 
does not exist a covariant curvature term of the mass dimension one. 



n-l 



V^V^GoVpGo ■ Go 

' :gpud p - gp P d v - g up dA5 D , (A.3) 



D + 2 n (2ire) n \D 



V^GoV^GoVpGo ■ Go™ 2 

— — f V9 p + g w d v + g vp d p )b D . (A.4) 



D(D + 2)n(n- 1) (27re 
Here, we have used the result. 

V M V,V P G ■ Go" 

which is obtained, in the same way as ( |A.1| ), from the symmetry under the exchange of the 
indices of the singular part. 

Next, we proceed to the singular behavior of the products of the functions Go's and Gi's 
with four derivatives. They have the mass dimension four. So, the singular piece consists 
of the bi-scalar delta functions with the quadratic derivatives or those multiplied by the 
curvature tensors. 

We cite the results for the singular parts of the products of G 's and G%s. 

V^GoVpV.Go • Go"" 1 

1 (_i)n ^ f2 _ D 



(D + 2)(D + 4) n (2vre)™ I D 
- 4(^V P V CT + g p(T V^ v )5 D 

+ D(g w V u V a + g^V u V p + g up V p V a + g vtT V „V P )5 D ^ 



g^gpa^ 6 + —(gp P g VCJ + gpag„p)V 5 



1 (—1)™ /i ne 
+ 3(D + 2) n (2-KeY M R *»"> + R ^ 

+ {(jppRua + g^crRup + 9upRpa + gucrRfip)}S D 

1 (-!)>- 1) f/ 2e\ * 

6(L> + 2)(L> + 4) n (27re)4V D) 9ilv9pa 

+ ~^{9pp9va + 9pa9up)R - ^{9pvRpa + 9paRpv) 
+ D(9ppRpo + 9poRvp + 9vpRpu + 9vaRp,p)^5 D 

6L> (2vre 



+ Rn/w^^^ + 9 P aRpi)5 D , (A.6) 
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:9^9pa^ 2 S D - Ugpp9.ua + 9pa9v P )V 2 $ 



€ ( l) n f 2 „ a w2c-D 1/- - i - A \V72r-D 



(£> + 2)(D + 4) n(n - 1) (2vre) n ID' 

- {gpp^v^a + gpa^u^p + gup^tNa + <kxV M V p )5 

e (—1)™ A*" 6 

3£>(£> + 2) n(n - 1) (27re)" ^ R »<™ + ^) 

+ iSppRva + (JpaRvp + gupRpa + gyaRpp)}^ 

e (-1)" f 2 A A ~ 

+ 6( J D + 2)(D + 4) n (27re)4 D 9ftv9pa 

+ 7,{gppgva + gp a gv P )R + {gpvRpv + g P vR, 



D 
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2 \*J pp*JW 1 fJ,V IS /J / 1 \J fJ,IS /L/L» 1 ^7 fJU pV 

(.9 ppRya 9p<jRvp 9upRpcr 9v<jRpp) 

— l) n /i ne r 1 - 

~Xg pvRpa gpoRpv) 
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6n (2TTe) n {D 

— -^(gppR-ua + gpoRvp + gv P Rp a + guaRp P )^S D . (A. 7) 

These correspond to the overlapping divergences, when a part of the x derivatives of Go's is 
converted into those with respect to y. 

We also have to consider the products involving a function G\. The equations ( 2.18Q 



make it possible for us to reduce the derivative of G\ to the sum of Go's an d G\S with 
no derivative. Using them, we find the singular parts of the products including a function 
G\ and four derivatives to be proportional to the delta functions without any derivatives or 



curvature tensors. There are three types of products necessary for our calculation. 



V^V.GoVpV.d ■ Go"" 1 

(-1)" . . a d 

~ 2D (2neY g ^ 9pa 

1 (-1)" fi ne f 2 „ A A A , „ 

+ 2(D + 2) n (2vre)™l D ^9pa + g^gu<r + 9^9u P )d 
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(A.10) 



4p + 2) (27re) 2 V£>- 

So far, we have considered the singular behavior of the products of Go's or Gi's with only 
x derivatives. Now, let us extend the previous results to the products with x' derivatives as 
well as those with respect to x. 

We can convert the x derivative of the functions Go, G\ to that with respect to x', using 
an identity which uses the parallel displacement bi- vector ■ 



7m' M V m G, 



-G,V^ , * = 0,1. 



(A.11) 



As a result, the singular expansions (|A.1| )- (|A.5| ) and (|A.8|) - ([A.10|) are valid under the change: 

^GiVuf, i = 0,l. (A.12) 



V^Gj 



1p 1 v p' j 



On the other hand, we should pay attention to the singular parts proportional to the 
quadratic derivatives of a delta function. In ( |A.6|) and (|A.7|) , we need additional curva- 
ture terms as 



XGo^vV ■ l/V u G Q ^ a , ■ Go"" 1 

V M V p G V,V CT Go • Go"- 1 
1 (-1)™ /i ni 



D + 2 n (27re) n lD 
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Rya + 9uaRpp) + R, 
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(A.13) 



l/V»G V p/ • V,G ■ ■ G n - 2 

~ V^VpGoV^GoV^Go • G n ~ 2 

e ( — l) n M ne A D 

+ 2D(D + 2)n(n-l) (2vre)" (jR ^ p<t + ^ ~ ^ R ^) 5 ° ■ ( A - 14 ) 



Here, we have used the coincidence limits ( |2.22| ). 
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